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Diffusion length measurements were made In highly doped and radiation damaged 
III-V semiconductors using the technique of charge collection microscopy (sometimes 
known as electron beam Induced current (EBIC)). EBIC curves were plotted whITe 
using the SEM on a Tine scan mode. Values of the currents read from these curves 
were then equated to expressions obtained from the solution of the diffusion equa- 
tion for a thick sample. An extended generation function was used In order to 
account for the finite volume of the Induced minority carriers. The surface recom 
bination velocity was either treated as an unknown In a system of two Integral equa- 
tions, or measured directly using low accelerating potentials for the electron beam. 

With the emergence of III-V compounds In the field of photovoltalcs , It has 
become increasingly more important to have accurate methods for determining small 
(1 to 10 pm) diffusion lengths L. A reliable determination of the values of L 
Is necessary In device modeling, radiation damage studies, and device fabrication 
since It Is quite Important to be able to assess the damage to the electronic pro- 
perties of materials subject to certain processes. 


Figure 1 shows the configuration that has been used to make these measurements. 
The depicted geometry was chosen to permit direct measurements to be made on 
finished solar cells. The contacts used were ohmic ones made with evaporated thin 
films of gold. The current amplifier had fast response, low noise, low input Imped 
ance, and the gain was calibrated. The leads were shielded, and the circuit ground 
was separated from the SEM ground. 
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energetic electrons impinge on semiconducting material, electron hole 
created. The required Ionization energy is a function of the bandgap 
In GaAs it takes 4.5 eV of incoming radiation energy to create a minority 
If one assumes no surface recombination, and point generation of minority 
the collected current follows a simple exponential decay form: 
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where 


| 0 maximum current collected (PN junction) 

L diffusion length 

X distance from PN junction 

However, the surface recombination can be very large in IIl-V's, and Its effect 
(ref. 2) cannot be neglected In the measurement of L. Increasing the accelerating 
potential for the electron beam minimizes the effect of surface recombination 
velocity S. Unfortunately, this approach diminishes the resolution of this 
technique, since for small L's the electron range then becomes comparable to the 
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value of L. Figure 2 shows that at 30 kV - the most commonly used accelerating 
potential in routine SEM operation - the electron range R e is about 3 pm. 

For these reasons low accelerating potentials have been used to obtain the 
charge collection microscopy curves used in these measurements and calculations. 
The Integral solution of the two-dimensional diffusion equation with semi-infinite 
thickness and an extended generation function (ref. 3) in the form of a three- 
dimensional gaussian distribution have been used here. Figure 3 shows an experi- 
mental EBIC plot of the ratio of collected current to maximum current (in the 
depletion region). At any x 0 , this ratio can be expressed as (ref. 4) 



As can be seen one can measure or assume reasonable values for all the variables in 
equation (2) except for the surface recombination velocity S, an unknown whose 
effect is not negligible. Hence, two different approaches that allow the determina- 
tion of S were undertaken, so that equation (2) could be solved numerically, and 
the value of L could be extracted from the Integral form. 

Assume two different accelerating potentials (E-| and E 2 ) for the electron 
beam impinging on the same semiconductor. In functional form, the normalized cur- 
rent at a given x 0 can be rewritten as 
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Figure 4 shows two of the experimentally obtained EBIC curves at the different 
potentials. The same x 0 is used, so one can assume the same value for the dif- 
fusion length. This applies even in the case of graded doping or other nonuniform- 
ities. At the same x 0 , one can also assume that S will be the same, even for 
different values of the recombination across the surfaces of the samples analyzed. 
This allows the treatment of equations (3) and (4) as two Integral equations with 
two unknowns. 
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The value for S can then be obtained by using an iterative process, where an 
Initial value for S Is guessed. Holding S constant In equation (3), an L Is 
found that satisfies the condition 


cc(calculated) 


cc(measured) 


< TOLERANCE 


(5) 


This value for L Is then used In equation (4), where S Is next varied to 
satisfy equation (5). This process Is repeated as necessary until an L and an S 
are found which satisfy both equations (3) and (4). Numerical Integration was done 
using the Romberg method where the upper limit was Increased until the last two com- 
puted Integrals differed by a negligible value. The roots (values of L and S) were 
searched by using the Regula Falsi Method. The Integral form for the complementary 
error function was used. Figure 5 shows a flowchart for the numerical calculations 
performed here. 


The second method uses less computer time but requires a more sophisticated SEM. 
It makes use of the result (ref. 5) 
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which allows a more direct determination of the value of S while the sample Is Inside 
the SEM specimen chamber. In order to use equation (6) and obtain accurate values, one 
must have low voltage capabilities and the ability to vary the beam accelerating poten- 
tial without changing the total beam current. Figures 6 and 7 show the determination 
of S/D for the devices that were analyzed here. The GaAs diode was P on N. The 
N region was silicon doped, with a carrier concentration of about lxlO 18 . The junc- 
tion was very abrupt. The InP solar cell had a P-type base, doped at about 1x10'^. 

The cell had been subject to ltP^/cm? 10 MeV proton Irradiation. 


Figures 8 and 9 show the measured values for diffusion length, as a function of 
distance from the PN junction, for the same devices. The spread In the data points 
from the different accelerating voltages (which Ideally, would coincide for a given 
x 0 ) has been used to assign a value to the uncertainty. The reported value has been 
chosen as the L that Is reached asymptotically as x 0 gets farther from the junc- 
tion, since other workers (ref. 6) have observed that the measured value of L Is more 
reliable If a larger x 0 Is used. 

In summary, accurate evaluations of diffusion lengths for heavily to moderately 
doped III-v semiconductors and/or radiation damaged solar cells have been made possible 
by using the experimental and numerical techniques described. 
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FIGURE 5. - FLOWCHART FOR NUMERICAL CALCULATIONS. USER INPUTS X , NORMALIZED 
CURRENTS, AND ELECTRON BEAM POTENTIALS FOR TWO CURVES. ° 


191 









